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Light sterile neutrinos mixing with the active ones have been recently proposed to solve different 
anomalies observed in short-baseline oscillation experiments. These neutrinos can also be produced 
by oscillations of the active neutrinos in the early universe, leaving possible traces on different cos- 
mological observables. Here we perform an updated study of the neutrino kinetic equations in (3+1) 
and (2+1) oscillation schemes, dynamically evolving primordial asymmetries of active neutrinos and 
ft j' taking into account for the first time CP-violation effects. In the absence of neutrino asymmetries, 

eV-mass scale sterile neutrinos would be completely thermalized creating a tension with respect 
• to the CMB, LSS and BBN data. In the past literature, active neutrino asymmetries have been 

invoked as a way to inhibit the sterile neutrino production via the in-medium suppression of the 
sterile-active mixing angle. However, neutrino asymmetries also permit a resonant sterile neutrino 
production. We find that if the active species have equal asymmetries L, a value \L\ = 1CF 3 is 
required to start suppressing the resonant sterile production, roughly an order of magnitude larger 
than what previously expected. When active species have opposite asymmetries the sterile abun- 
O I- dance is further enhanced, requiring an even larger \L\ ~ 10~ 2 to start suppressing their production. 

In the latter case, CP-violation (naturally expected) further exacerbates the phenomenon. Some 
consequences for cosmological observables are briefly discussed: for example, it is likely that moder- 
ate suppressions of the sterile species production are associated with significant spectral distortions 
of the active neutrino species, with potentially interesting phenomenological consequences especially 
for BBN. 
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O , I. INTRODUCTION 

In recent years a renewed attention has been paid to light (m ~ 0(1) eV) sterile neutrinos mixing with the active 
r"-; ones (see Q for a recent review). In particular, sterile neutrinos have been proposed to solve different anomalies 
observed in short-baseline neutrino experiments, notably in the — > T> e oscillations in LSND Q and MiniBoone H 
experiments, and in T> e and v e disappearance revealed by the Reactor Anomaly [J] and the Gallium Anomaly [5j, 
respectively. Scenarios with one (dubbed "3+1") or two ("3+2") sub-eV sterile neutrinos [6l-tl0| have been proposed 
to fit the different data. 

Cosmology provides an important arena to test these scenarios. In fact, neutrinos are abundantly produced via weak 
interactions in the hot (temperature T ^ 1 MeV) primordial cosmic soup. The mostly sterile mass eigenstate(s) can 
be produced via oscillations and modify cosmological observables fllTflij . If these additional states are produced well 
before (active) neutrino collisional decoupling, they acquire quasi-thermal distributions and behave as extra degrees of 
freedom at the time of big bang nucleosynthesis (BBN). This would anticipate weak interaction decoupling and leading 
to a larger neutron-to-proton ratio, eventually resulting into a larger 4 He fraction. The non-electromagnetic cosmic 
radiation content is usually expressed in terms of the effective numbers of thermally excited neutrin o sp ecies N e f{ . The 
Standard Model (plus active neutrino oscillations) expectation for this parameter is N c g = 3.046 [15|, a result which 
is only marginally modified even accounting for new interactions between neutrinos and electrons parameterized by 
four-fermion operators of dimension six allowed by laboratory constraints [l6| ■ If the additional degrees of freedom 
are still relativistic at the time of cosmic microwave background (CMB) formation, the same parameter N g can 
be constrained by a detailed study of the CMB angular power spectrum, especially when combined with other 
cosmological probes. Quite an excitement (see e.g. [13, [HI) nas been stimulated by the result in the current best 
fit of WMAP, SDSS II-Baryon Acoustic Oscillations and Hubble Space Telescope data, yielding a 68% CL range on 
N c fi = 4.34lg gg [lj| in the assumption of a ACDM universe. Once accounting for the parameter degeneracy in the 
determination of the first angular peak properties, which can be adjusted by choosing different combinations of other 
cosmological parameters, it turns out that such results are almost completely due to the large-£, damping tail of the 
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CMB spectrum, as even more clear when adding ACBAR [2(| and ACT (2l[ small scale data (see [22| for a pedagogical 
account). The reliability of BBN constraints is plagued by systematics in the determination of primordial abundances. 
Yet, within conservative but reasonable assumptions, standard BBN calculations do not allow for N e s larger than 
about 4.1 at 95% CL [23[ with only a weak, statistically non-significant preference for a larger-than-standard value 
of iV e ff. In order to allow for sufficiently large effects at CMB epoch while accommodating BBN constraints, some 
authors have envisaged for example the introduction of relatively large neutrino-antineutrino asymmetries. This in 
order to (partially) compensate for the effect of N e g on He by a counter-effect due to v e — V e distributions in weak 
rates [241 ]. On the other hand, the CMB preference for a large N e g usually comes with a price: for example, a tension 
with cluster determination of dark matter abundance has been noticed in (22j. Also, the basic tenet that these mostly 
sterile i/'s behave essentially as radiation at CMB epoch is untenable: laboratory data require one or more relatively 
massive (m ~ 1 eV) extra states. CMB data in combination with LSS ones are a particularly sensitive probe of 
massive neutrinos (for a review see e.g. [25|,[26|]). At face value, when accounting for the fact that these extra species 
are massive, the scenario hinted at by laboratory data is actually disfavored by cosmology, unless rather radical and 
contrived cosmological model modifications are introduced [HI, H3, l27l - [29l | . 

Given this partially contradictory situation and the existing laboratory anomalies, it is of paramount importance 
to study the physical conditions under which the sterile neutrino production actually takes place, as preliminary step 
to any phenomenological consideration. As already mentioned, sterile neutrinos are produced in the early universe 
by the mixing with the active species. Therefore, in order to determine their abundance it is necessary to solve 
the quantum kinetic equations describing the active-sterile oscillations (30L l3lj ] . This problem has been studied in a 
long series of papers (see e.g. (32| - [53j ). finding a broad range of possible outputs depending on the sterile neutrino 
masses and mixing angles with the active species. Since the solution of the non-linear neutrino evolution equations 
is numerically challenging, different approximations have been adopted. In particular, most of the previous studies 
solve the equations in a simplified (1+1) scenario in which only a (mostly) active neutrino mixes with a (mostly) 
sterile one. Recently, also multi-flavor cases have been presented [48|,[52[. In particular, active-sterile oscillations in a 
(3+2) scheme have been studied 52]. It has been found that, for the mass and mixing parameters needed to describe 
the short-baseline anomalies, in the standard scenario sterile neutrinos would be completely thermalized in the early 
universe, creating an unwelcome tension with cosmological observations as mentioned above. 

On the other hand, a possible escape-route to reconcile sterile neutrinos with cosmological data consists in the 
inclusion of a primordial asymmetry between neutrinos and antineutrinos (4lj 

L = — - Ue . (1) 

rij 

In principle, one would expect the lepton asymmetry to be of the same order of magnitude of the baryonic one, 
V = ( n B — n B)l n i — 6 x 10~ 10 . Indeed, to respect the charge neutrality the asymmetry in the charged leptons must 
match the above number to a high degree. However, since neutrinos are neutral the constrains on L are quite loose, 
allowing also \L\ ~ 1CP 2 — 10 _1 [5J-[6fJ. Moreover, there are models for producing large L and small r\ [6ll-t63l]. 

A neutrino asymmetry implies an additional "matter term potential" in the equations of motion. If sufficiently large, 
one expects this term to block the active-sterile flavor conversions via the in-medium suppression of the mixing angle. 
However, this term can also generate Mikheev-Smirnov-Wolfenstein 64] (MSW)-likc resonant flavor conversions among 
active and sterile neutrinos. In particular, in a recent (1+1) study [531 ] it has been found that for sterile neutrinos 
with parameters preferred by the laboratory hints, a neutrino asymmetry L = 10 -2 would strongly suppress their 
production. This would reconcile them with the cosmological observations. In [50| the question was addressed of 
how large should be the value of L in order to have a significant reduction of the sterile neutrino abundance. The 
authors solved the equations of motion in a simplified (3+1) scheme inspired by LSND, finding that L ~ 1CP 4 was 
enough to relieve the tension between sterile neutrinos and cosmology. However, this result has to be taken cum 
grano salis. Indeed, the authors fixed the lepton asymmetry as an initial condition taken constant during the flavor 
evolution. Nevertheless, this quantity is expected to dynamically evolve due to the flavor conversions. Moreover, they 
solved the coupled equations of motion by effectively reducing the degrees of freedom via the constraint of neglecting 
resonant transitions between sterile and active neutrino species in the antineutrino sector, that would be possible for 
the negative neutrino asymmetries they considered. 

only for neutrinos, neglecting the antineutrino sector, in which resonant conversions with active neutrinos would 
occur for the negative lepton asymmetries they considered. 

The purpose of our work is to revisit the thermalization of sterile neutrinos in the early universe in the presence 
of primordial neutrino asymmetries, taking a complementary approach to most of other studies. In particular, many 
investigations have focused on large scans of the sterile neutrinos parameter space, often neglecting some physical 
effects in order to keep the problem computationally manageable. Here, we rather stick to a benchmark set of best fit 
value parameters for the sterile state "suggested" by laboratory experiments, but we go beyond most approximations 
used in the previous studies. In particular, we shall consider (3+1) and (2+1) schemes inspired by the recent fits of all 
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the short-baseline, reactor, and solar neutrino anomalies. The inclusion of more than an active neutrino that mixes 
with the sterile state allows to explore effects that are not possible in a simplified (1+1) scenario. For example, having 
more than one active neutrino one can study both cases in which the neutrino asymmetries are equal or opposite 
between the active species. Moreover, one can take into account more than one mixing angle between the active and 
the sterile neutrinos. CP violating effects in oscillations thus become a natural possibility, which we also consider here, 
to the best of our knowledge for the first time in this context 1 . We decide to devote our work to a detailed study of 
all these effects still unexplored. We find that each of them could have a relevant impact in the determination of the 
final abundance of sterile neutrinos. Generically, we realize that "masking" cosmological consequences of additional 
sterile neutrinos is harder than previously deduced in simplified treatments. Here we do not aim at deriving detailed 
cosmological predictions, rather the correct (and reliable) qualitative physical trend. In this spirit, we can content 
ourselves with a momentum-averaged description, but otherwise we deal with the complete problem, i.e. we solve 
essentially the exact equations of motion as opposed to simplified ones studied in the past literature. 

The plan of our paper is as follows. In Sec. II we introduce the (3+1) neutrino mixing framework we will use as 
benchmark for our study. In Sec. Ill we describe the active-sterile neutrino flavor evolution in the early universe. In 
particular, we introduce the equations of motion for the neutrino ensemble. Then, we present the average-momentum 
approximation we use to solve numerically the neutrino evolution. Finally, we compare the different interaction 
strengths that enter the equations of motion. In Sec. IV we present the results of the sterile neutrino production in 
the (3+1) scenario for different values of the primordial neutrino asymmetries, taken equal among the different active 
flavors. In Sec. V we calculate the sterile neutrino abundance in different cases in a (2+1) scenario, where the active 
sector (v e ,v^) is associated with the atmospheric mass-square splitting Am^ tm and with the 1-3 mixing angle #13. 
At first, we take into account both the active-sterile mixing angles 9 es and 6^ s . In this situation we describe both 
cases with equal and opposite initial neutrino asymmetries among the active species. We also consider the effects of 
CP violation in the sterile sector. Then, we present a case in which only electron neutrinos mix with the sterile ones. 
In Sec. VI we discuss the impact of the sterile neutrino production on the neutron/proton ratio (n/p) in the early 
universe. Finally, in Sec. VII we comment about future developments of our study and we conclude. 



II. (3+1) NEUTRINO FRAMEWORK 

Short-baseline neutrino oscillation data have been fitted adding to the usual three active neutrino species one (3+1 
mixing scheme) or two (3+2 mixing scheme) massive sterile states [H-IH . In spite of the presence of a tension in the 
interpretation of the data the 3+1 neutrino mixing scenario is attractive for its simplicity. It is also less likely to 
lead to a strong exclusion by cosmological arguments. Therefore, in our work we will consider only this extension of 
the three-neutrino mixing framework. In such a four-neutrino mixing scheme, the flavor neutrino basis is composed 
by the three active neutrinos v^v^, v T and by a sterile neutrino v s . The flavor eigenstates v a are related to the mass 
eigenstates i/j (i = 1, . . .4, ordered by growing mass) via a unitary matrix U through p], [66| 

v a = W m v, , Ulfi =U^U = I , (2) 

By neglecting for the moment the presence of arbitrary phases, also responsible for CP violation effects, the matrix 
U can be parameterized as a product of 4x4 Euler rotation matrices -Ry acting in the (i, j) mass eigenstate subspace, 
each characterized by a mixing angle Oij. Following e.g. [67| one can write 

U = -R34-R24-R23-R14-R13R12 , (3) 

where we order the flavor eigenstates in such a way that if all angles are vanishing we have the correspondence 
(i> e , Vp, v T , u s ) = (y\, vii V3, va)- In the limit where the three mixing angles On vanish, the above matrix reduces to 

lim U =( U ^i lM °A , (4) 
i4 ^o \ 1 J 

where U is the 3x3 unitary mixing matrix among the active neutrinos defined in terms of three rotation angles 9^ , 
ordered as for the quark mixing matrix 68]. In the following we fix the values of these three mixing angles to the 
current best-fit from global analysis of the different active neutrino oscillation data [69| (see also [7(|), i.e. 

sin 2 6»i2 = 0.307 , 



1 In the absence of sterile neutrinos, the impact of CP violation onto cosmological active neutrino asymmetries has been studied in |65H . 
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sin 2 e»23 = 0.398 , 

sin 2 6»i3 = 0.0245 . (5) 

We remind that the early hints for a "large" value of 6*13, suggested by the long- baseline v^-v e experiments [7ll. W\ 
and Double Chooz reactor experiment [73| in combination with the global analysis of the neutrino data [lH, 175} . have 
been recently confirmed by the measurement of the Daya Bay [76| and Reno [73 reactor experiments. We neglect CP 
violating effects in the active sector. The corrections to the 3x3 "active" neutrino sub-matrix of U is only second 
order in the mixing angles of the sterile state. We shall assume (as also done in phenomenological studies) that at 
most two mixings of the fourth neutrino to the three active ones are non vanishing, namely we put IA T ^ = 0. We take 
from fits of the short-baseline data the values @ 



sm6 es ~ |W e4 | = V0.025 , 

sin^ s ~ |w m4 | = Vomz , (6) 

where the equality holds (in our parameterization) up to corrections quadratic in the mixing angles of the fourth state. 
It is worth noting that the "sterile" mixing angles result of the same order of 6*13 . Any quantitative interpretation of 
"anomalies" in terms of mixing with steriles should thus be made to account for oscillations among the active states. 
The Av mass spectrum is parameterized as [78j 

M 2 = diag(m 2 , to 2 , m§, m\) = diag (-^Am* ol , + ^Ato 2 o1 , Am 2 tm , Am^ t J , (7) 

where the solar and the atmospheric mass-square differences are given by (69j 

Am 2 ol /eV 2 = 7.54 x 10~ 5 , 

Am 2 tm /eV 2 = 2.43 x 10" 3 , (8) 

respectively. Note that here and throughout we assume normal mass hierarchy, i.e. Am 2 tm > 0. The sterile-active 
mass splitting from the short-baseline fit in the 3+1 model is given by Q 

Am 2 t /eV 2 = 0.89 . (9) 

Therefore, it results in a clear hierarchy among the mass differences, i.e. Am 2 t 3> Am 2 tm 3> Am 2 ol . When restricting 
to the (2+1) cases, no 4 x 4 formalism is needed. The conventions are formally equivalent to the familiar three active 
neutrino mixing angles, as well as the standard parameterization of the Dirac CP phase [68j], with only the numerical 
values for the mixing and mass splittings to be changed (see Sec. [V]for details). 

III. NEUTRINO FLAVOR EVOLUTION IN THE EARLY UNIVERSE 

A. Equations of Motion 

Following [EiJ], in order to describe the time evolution of v — V ensemble in the early universe, it proves useful to 
define the following dimensionless variables which replace time, momentum and photon temperature, respectively 

x = ma y = pa z = T^a, (10) 

where to is an arbitrary mass scale which can be put e.g. equal to 1 MeV. Note that the function a is normalized, 
without loss of generality, so that a[t) — > 1/T at large temperatures, T being the common temperature of the particles 
in equilibrium far from any entropy-release process. With this choice, a -1 can be identified with the initial temperature 
of thermal, active neutrinos. 

In order to characterize the active-sterile neutrino oscillations we describe the neutrino (antineutrino) ensemble in 
terms of 4 x 4 density matrices g (g) 2 
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(11) 



2 Note that in natural units q and g are dimensionless variables. 
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In terms of g and g the Equations of Motion (EoMs) for the neutrino ensemble assume the form [3CJ, |3l|, [5 
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In the previous expressions H denotes the properly normalized Hubble parameter, namely 

,•2 



m m 



1 8w e(x, z(x)) ( m \ j 8tte(x, z(x)) 
= \~Mp~i) ' 



3M 2 Pl 



(12) 

(13) 
(14) 

(15) 



where the total energy density and pressure of the plasma, e and P, enter through their "comoving transformed" 
£ = e{x/m) i and V = P(x/m) 4: respectively. Since for most of the temperatures we are interested in, electrons and 
positrons are the only charged leptons populating the plasma in large numbers, to a very good approximation the 
total energy density can be expressed as the sum 



e(x, z(x)) 



where 



— z (x) . 
15 v ; 



1 



dyy 3 [f F D{y/z{x) - 4> e ) + f FD (y/z(x) + <j> e )] 



7t? 
60 



z 4 (x) , 



2^2 J d Vy 3rTr [s(x,y) + g{x,y)] = ^jc- N e$- 



(16) 

(17) 
(18) 
(19) 



Note that due to the range of temperature T considered we have safely assumed massless that, due to the fast 
electromagnetic interactions, have a Fermi-Dirac distribution f F o{y I 'z{x)^<j> e ) = l/(exp(y/z(a;)=F0 e ) + l) respectively. 
The reduced electron chemical potential cj) e is in principle a dynamical variable that requires a further equation (the 
electric charge conservation) in order to be evolved consistently. However, for our purpose electrons are only important 
when their energy density is dominated by pairs, rather than by the e~ excess due to the baryon asymmetry, and 4> e 
can be put equal to zero. 

The first term on the r.h.s. of the EoMs (fT2"j) and (|13l) is responsible for the vacuum neutrino oscillations. In the 
second term, the diagonal matrix Ef related to the energy density of charged leptons under the previous assumptions 
takes the form 



diag(e e , 0, 0, 0) = diag 



7t? 
60 



z 4 (x), 0,0,0 



Moreover we have 



N, 



1 

1 

2^2 



dyy 2 {G s (g(x,y) - g(x,y))G s + G s Tr [(g(x, y) - g(x, y))G s ]} , 
dyy 3 G s (g(x,y) + g(x,y))G s . 



(20) 

(21) 
(22) 



These terms make the EoMs non-linear and are the main numerical challenge in dealing with this physical system. 
Note that the matrix N„ is related to the difference of the density matrices of neutrinos and antineutrinos, while E„ 
is related to their sum. The matrix G s = diag(l, 1, 1,0) in flavor space contains the dimensionless coupling constants. 
We remark that in the presence of more than one active species, the matrix also contains off-diagonal terms. The 



6 



last term at r.h.s. of Eqs. (IT21) and (IT3|) is the collisional term proportional to G 2 F . We will present an approximate 
expression of this term in the Sec. IIII Bl Finally, Eq. ((TJJ basically provides an evolution equation for the quantity 
z(x). On the other hand, even when a fourth neutrino is populated in the plasma via early oscillations (before the 
active neutrino decoupling), the correction with respect to the initial value of e is at most of the order of ~ 10%. 
Consistently with other papers in this field we neglect this small effect in this exploratory study, although we note 
that it should be accounted for in more accurate predictions of cosmological observables. This also implies that we 
can keep z as constant and equal to 1, which is a good approximation in the epoch considered in this work. As a 
consequence, e is not dynamical and assumes a numerical value almost equal to 3.54. A final comment is that we omit 
the familiar refractive term due to the matter asymmetry (i.e. oc (n e - — n e +) cx rf) from the EoMs. Albeit, it may, 
in principle induce resonances, for the parameters of interest here the resonance would fall at such high temperature 
that the system is still collisionally dominated and no coherent process is effectively taking place. 



B. "Average momentum" approximation 



In the presence of continuous neutrino momentum distributions, to solve the full set of EoMs (|12l) and (|13l) turns 
out to be a computationally demanding task. In order to perform a more treatable numerical study of the flavor 
evolution for different neutrino asymmetries, which is able to catch the main features of the more involved complete 
computation, in the following we will restrict ourselves to a average-momentum approximation, based on the Ansatz 
(and similarly for antineutrinos) 



Q{x,v) -> fFD(y) p{x) . 



(23) 



By this assumption, and in absence of asymmetries, at equilibrium one would simply get p = I. In terms of Eq. ([23 
the set of EoMs (|T2"j) and (flU)) can be rewritten as 
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(24) 



(25) 



where by definition (g(y)) = J °° y 2 g(y) /fd(j/) dy/ J °° y 2 /f.d(2/) dy. According to this notation (y) = 3.15 and 
(1/y) — 0.456 7^ l/(y)- The non-linear terms N„ and of Eqs. (|21I22[) assume the form 



3C(3) 
4tt 2 

7 7T 2 



{G s (p-p)G s +G s Tv[(p-p)G s }}. 



E " = 8 30 Gs( ' P + P ' )Gs ' 



(26) 
(27) 



Note that only the "active" 3x3 submatrix of the whole density matrix enters the above interaction terms. Moreover, 
by using the approximate form [5(| for the collisional terms in Eqs (fT2"|) and (|13[) . we get the expressions 



C[ P ] = — 2 G 2 F ; 
C[p] - -\G 2 F >, 



'({S 2 , p - 1} - 2S( P - l)S + {A 2 , [p - I)} + 2A(p - l)A) , 
'({S 2 , p - 1} - 2S{p - l)S + {A 2 , (p - I)} + 2k(p l)A) . 



(28) 



The above expressions have been derived assuming null neutrino lepton asymmetries. However, since in our study 
we will restrict ourselves to \L\ < 10~ 2 , the correction induced would be negligible (see, e.g., [43|]). In flavor space, 
the matrices S,A write S = diag(<?f, g T sl 0) and A = diag(g^, g%, g T a , 0), respectively, and contain the numerical 
coefficients for the scattering and annihilation processes of the different flavors. Numerically one finds [33] 



{9lf 



3.06 
2.22 



0.50 , 

0.28 



(29) 
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The initial conditions for the density matrix p are written 

/ 4 4 4 

diagf 1 + -L e ,l + -£ M ,l + -i T ,0 



An = diag fl-|i e ,l -1l m> 1-|l t ,0^ , (30) 

where the neutrino asymmetries in the different ffavors are related to dimensionless chemical potentials = p v a 
through 

12C(3) V s " ' n 5 

Note that the expression in Eq. (IBTfl) is only valid at leading order in L. 

It is convenient to have an estimate of the different dimensionless factors multiplying p and p on the r.h.s. of 
Eqs. (|24[) and (|25[) . respectively. The vacuum oscillation term is proportional, apart from a matrix whose coefficients 
are 0(1), to the quantity 

1 \ „ „ . „ 13 



^ = 7^tU, + ^t). (31) 



f2 vac = — = ( - ) = 2.3 x 10~ iJ 5- = , (32) 

2ffm 2 W V eV 2 J H V ' 



Taking into account the e + e pairs only, the matter potential in Eqs. (|24|) and ([25]) except for the different sign for 
neutrinos and antineutrinos, is proportional to 

o matt = g^Mg^!g f M = 2 . 4 x 10 -*> i . (33) 

3™^ 60 Va; 4 ^/ a: 4 ^ 
The neutrino-neutrino interaction strength gives two terms, respectively proportional to 

0^ = ^^^ = 1.5x10- 1 



8V2(y)G F m 4 7tt 2 / 1 \ 1 /m ff \ 2 cos 2 <V 



Finally, the collisional term is proportional to 
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G F TO _ R o v 1n -23 1 



^, = ^=6.8X10 --j,. (35) 

In order to get an idea of the strength of the different interaction terms, in Figure 1 we plot as a function of 
the temperature T, Q vac (solid curve), f2 ma tt (long-dotted curve), i7 asy x A e (dashed curve), ri S ym x S ee (short- 
dotted curve), n co ii x [(gf) 2 + (s^) 2 ] (dash-dotted curve). Here we use as mass square difference Am 2 t , A e = 
2(/5 e e — Pee) + (p^[i — /W) + (/3rr — p TT ) = 8/3(2£ e + + L T ), where for illustration we fixed L e = = L T = 10~ 4 . 
Finally E ee = (p ee + p ee ) ~ 2. 

From the Figure above, one realizes that the system remains collisional down to a few MeV, when the collision over 
Hubble rate drops below 1 on the ordinates. The collisional term also dominates over the vacuum oscillation term at 
T > 20 MeV, thus breaking the coherence between different neutrino flavors and preventing significant oscillations. 
The refractive terms can induce MSW-like resonances between the actives (a = e, p, r) and sterile state when, in 
the limit of only one mixing angle between the active and the sterile neutrinos, one of the following conditions is 
satisfied [H 

^vac COS 20 as — £7 aS y A a -\- ^sym^aa ~i~ ^mat fc*l V , 

fi vac COS 26 as + Oasy A + O sym S ao + O mat = for 77, (36) 

where a — e, p,r and the definitions of A M , A r and E Tr are respectively of the same form to the ones of A e 

and E ee given before. From these equations, we obtain that in absence of lepton asymmetries (A a = 0) the resonance 
condition cannot be satisfied neither for the i/'s nor for 77's, given the hypothesis that the sterile state is heavier than 
the active ones. Instead, when 57 asy A a is the dominant term, as in the cases we will consider in the following, resonance 
conditions can occur for A a > in the v sector and for A a < in the V one. In particular, in Fig. [1] the resonance 
occurs around T ~ 3 MeV. We will also show that, as a consequence of the dynamical nature of the asymmetries, 
A a can rapidly change sign so that both sterile neutrinos and antineutrinos get populated. This phenomenon is thus 
qualitatively different with respect to the familiar MSW resonant conversion. Resonances can also take place in the 
active sector at lower temperatures. However, since active neutrino distributions are expected not to depart too much 
from their equilibrium values, their effect is sub-leading. 
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FIG. 1: Neutrino collisional and refractive rates (normalized in terms of the Hubble rate) vs. temperature T. In particular, we 



show f^vac (solid curve), f2 m att (long-dotted curve), 
(#s ) 2 ] (dash-dotted curve). Here we use Am;? 
quantities). 



(dashed curve), fi sym X S ee (short-dotted curve), f2 C oil x [(ffs) + 
32L/3 with L = 10~ 4 and E ee = 2 (see the text for more details on these 



IV. (3+1) RESULTS 

In order to calculate the sterile neutrino abundance in the 3+1 scenario, described in Sec. II, we numerically 
solved the EoMs [Eq. ([25]) ]. using a Runge-Kutta method for the equations written in the variable x = m/T and 
evolved in the range x £ [10 -2 , 1.0]. We take 10 5 steps in log(x) in the integration interval. We consider initial 
neutrino asymmetries L = L e , = = L T < 0. We checked that the results presented in the following do not change 
considering positive asymmetries. In Fig. [2] we show the evolution of the diagonal component of the density matrix 
p ss for sterile neutrinos in function of the temperature T for different initial lepton asymmetries, namely L = 
(solid curve), L = — 10~ 4 (dashed curve), L = — 10~ 3 (dotted curve), and L = — 10~ 2 (dash-dotted curve). As 
expected from the previous literature, in absence of lepton asymmetries sterile neutrinos are copiously produced at 
T < 30 MeV until they reach p ss = 1. Instead, including a non-zero initial lepton asymmetry the effect is to suppress 
the sterile neutrino production as long as |^asy| 3* |^vac|- However, these two functions have opposite dependence 
on the temperature and at some time they will cross. Sterile neutrinos are then produced "resonantly" , albeit with 
a non-linear, dynamical resonance condition which is itself influenced by the evolution of the system. Increasing the 
lepton number asymmetry the position of the resonance moves towards lower temperatures, where the resonance is 
less adiabatic. Indeed, the adiabaticity parameter scales as ~ T, as shown in 79]. As a consequence, the sterile 
production is less efficient increasing \L\, as results from Fig. [5] In particular, asymmetries greater than \L\ = 10~ 3 
are required in order to achieve a significant suppression of the sterile neutrino production. Also, the asymmetric 
term changes sign and thus the resonance can take place in both neutrino and antineutrino sectors, which turn out 
to be populated almost equally. 

Our result implies that in order to suppress the sterile neutrino production one needs a lepton asymmetry greater at 
least by an order of magnitude with respect to what found in a previous study on the subject [501 ] . This discrepancy is 
due to the fact that in their work the authors followed the flavor evolution only for the neutrinos, choosing a negative 
value of lepton asymmetry kept constant. In this way, they missed resonant effects that would have occurred in 
the antineutrino sector. Then, the lepton number can only suppress the flavor evolution. Therefore, in their study 
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FIG. 2: (3+1) scenario. Evolution of the density matrix element p ss in function of the temperature T. We consider L = L e = 
£p = L T . The solid curve corresponds to L = 0, the dashed curve to L — — 10~ 4 , the dotted curve to L = — 10 -3 and the 
dash-dotted one to L — — 10 -2 . 



L = — 1CP 4 was enough to block the sterile neutrino production. 

Caution should also be taken when interpreting the results shown for p ss into an effective increase of the neutrino 
degrees of freedom in the early universe, usually parameterized via N^g. In fact, according to the definition reported 
in Eq. (|19[) . the latter variable is sensitive to the trace of the neutrino plus antineutrino density matrix. A late 
conversion of some active state into a sterile one after the neutrinos have undergone collisional decoupling is in fact 
conserving the overall number of neutrinos (albeit some cosmological consequences, such as those for BBN, may be 
typically more dramatic, as briefly discussed in Sec. IVI[) . This is shown in Fig. [3J reporting the evolution of N e s for 
the cases corresponding to Fig. [5] Note that for no or small asymmetry, for the parameters chosen the active-sterile 
oscillations take place early enough that the depleted active states are rapidly repopulated collisionally. Thus N e g 
effectively increases to 4. On the other hand, for \L\ = 10~ 3 the conversion takes place around the decoupling time, 
and the repopulation is only partial, with a difference between N e g and p ss of about 0.1 units (compare Fig. [5] with 
Fig. [3]). Finally, for \L\ — 10~ 2 only a negligible fraction of the converted active neutrinos are repopulated, despite 
the fact that about 10% of a "thermal-equivalent" sterile state has been produced. Since of large asymmetries the 
temperature at which production starts depends on when the equality f2 vac = f2 asv x A e takes place, there is a quite 
strong dependence of the signatures from the exact values of the active neutrino mass, mixing, and the initial value 
of ILL 



V. (2+1) RESULTS 



In our study we consider initial distributions for active neutrinos close to their equilibrium ones. Therefore, the 
oscillations among the three active species have a sub-leading role for the evolution of the sterile neutrinos. At this 
regard, we calculated the flavor conversions in the same cases as before, considering (2+1) sub-sectors with the active 
mixing associated with (Am^, On) and (Am atm , #13), respectively. For the cases we compared, we find results very 
similar to the ones presented in the previous section. Therefore, in order to speed-up the numerical calculations we 
decide to continue our explorations of sterile neutrino production in different cases, referring to (2 + 1) scenarios 
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FIG. 3: (3+1) scenario. Evolution of the effective number of degrees of freedom A^fr for the cases corresponding to Fig. 



associated with (Am^ tm ,0i 3 ) active sector. 



A. L e = L M , ipcp = 

In the following we consider different (2+1) cases with non-zero 6 es and Q^s given by Eq. (j6|). In the left-upper 
panel of Fig. @] we represent the case with L = L e = L^. The solid curve corresponds to L = 0, the dashed curve to 
L = — 10~ 4 , the dotted curve to L = — 10 -3 and the dash-dotted one to L = — 10~ 2 . This case is manifestly close 
to the (3+1) scenario shown in Fig. [2] In order to clarify the dynamics of the sterile neutrino production, in the left 
panels of Fig. [5] we plot in function of the temperature, the evolution of the neutrino asymmetries Ap a = p aa — p aa 
for the v e (solid curve), (dotted curve) and v s (dashed curve). Since Ap a typically presents very fast oscillations, 
for the sake of the clarity we plot its value averaged over ten steps in T . In the right panels we show the evolution of 
the vacuum term fl vac (solid curve) and of the f2 a sy x A e term (dashed curve) for the same cases of the left panels. 
The crossing between these two curves at non-zero L determines the position of a e-s resonance. 

Starting with the case L = 0, we see that L e = —2L^ = —2L S ~ few x 10~ 5 can be dynamically generated at the 
onset of the flavor conversions (at T < 80 MeV). Since the active asymmetries are opposite, they tend to decrease 
reaching flavor equilibrium (L = 0) at T ~ 10 MeV. At T < 30 MeV, when collisional rates slow down enough (see 
Fig- [U i sterile neutrinos are produced without any hindrance (see Fig. 2]). 

We pass now to the cases with non-zero initial neutrino asymmetries. In this situation, since 9 es and 9^ s are non- 
vanishing, both the active states can have resonances with the sterile one. Moreover, since for our choice 6 es ~ 6^ s , 
the evolution of Ap e and Ap^ is very similar. 

In the case with initial L = — 10 -4 the production of v s starts at T ~ 10 MeV (Fig. 2]) when an active-sterile 
resonance occurs. Also in the other two cases with L = — 10~ 3 and L = — 10~ 2 the position of the resonance coincides 
with the onset in the rise of p ss in Fig. |U However, as commented before, the lower the resonance temperature, the 
less adiabatic the resonance. Therefore, the sterile neutrino production is further inhibited. 
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FIG. 4: (2+1) scenario. Evolution in function of the temperature T of p sa for different initial neutrino asymmetries. Upper 
panels correspond to L = L e = L M , lower panels correspond to L = L e — ~ L M . The solid curves correspond to L = 0, the 
dashed curves to L — — 10 -4 , the dotted to L = — 10 -3 and the dash-dotted to L = — 10 -2 . Left panels show cases with no CP 
violation in the sterile neutrino sector, while right panels refers to ificp = it/2. 

B. L e = L M , )3CP = tt/2 

Fits for laboratory anomalies have been proposed which include CP violation effects in the sterile sector (see 
e.g. [80j|). Perhaps more importantly, whenever three or more neutrinos mix, CP- violating "Dirac phases" entering 
oscillations are naturally present in the theory. For both reasons, we find worthwhile to investigate the impact of 
CP-violation in our framework. 

For this purpose, we include an extra phase in the sterile- active mixing matrix [Eq. ([¥])]. formally in the same way 
the Dirac phase is introduced in 3 x 3 active neutrino mixing formalism. The inclusion of CP violating effects in 
the sterile sector could be potentially interesting, since it would generate an asymmetry among sterile neutrinos and 
antineutrinos. This could be transferred by oscillations into the active sector, having a feedback on the further growth 
of the sterile neutrino abundance. For defmiteness we consider tpcp — t/2. Note also that in the full (3+1) scenario — 
not to speak of the (3+2) scenarios with 2 sterile states — the number of CP-violating phases grows. Therefore, the 
present investigation is expected to be conservative in some respect. We first refer to the case with initial equal 
neutrino asymmetries among active species: L = L e = L^. The evolution of p ss is shown in the right- upper panel of 
Fig. |31 whose comparison with the CP conserving case shows that the suppression of the sterile neutrino abundance 
due to ifcp is sub-leading. Indeed, from Fig. [6] one sees that the growth of the dynamical neutrino asymmetries Ap a 
for different intial L, even if it is more irregular than in the case with ipcp — (Fig.[S]), it is qualitatively similar. This 
implies that this effect docs not significantly alter the flavor evolution. On the other hand, it is interesting to note 
that even in absence of an initial neutrino asymmetry the CP-violating mixing can create a "dynamical" asymmetry 
at relatively late times (down to decoupling temperatures) which is of the order to 10 -5 for the parameters used. 
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C. L e = — L M , ipcp — 

We pass now to consider the case in which the initial neutrino asymmetries in the active sector are opposite for 
v e and z/„, i.e. L = L e = —L„. In the absence of CP violation, this case is represented in the bottom-left panel of 
Fig. 2J For a non-vanishing initial L the sterile neutrino production is enhanced with respect to the previous case 
with equal asymmetries among the different flavors. Indeed, to achieve a significant suppression of the sterile species 
one needs an initial \L\ ~ 10~ 2 , i.e. roughly one order of magnitude larger than in the previous case. This behavior 
can be clarified looking at the evolution of the dynamical asymmetries Ap a shown for the different initial L in Fig. [7] 
We remark that since Ap e and Ap M have opposite sign, resonances can occur simultaneously in the neutrino and 
antincutrino channels. When these happen, they tend to produce flavor equilibrium between v e and v^. This leads to 
a vanishing final lepton number. When the neutrino asymmetry is destroyed, the sterile neutrinos can be produced 
without any hindrance. This explains the enhancement in the final p ss found in this case. 

D. L e = -Lp, ip C p = 7r/2 

We now consider the case with opposite initial neutrino asymmetries and (pep — 7r/2. The evolution of p ss in this 
case is shown in the bottom-right panel of Fig. 01 The production of sterile neutrinos is significantly enhanced with 
respect to the previous cases. In particular, also for an initial L = — 10~ 2 , the final abundance of sterile neutrinos is 
relevant. From Fig. [5] one sees that the flavor equilibrium between the electron and muon species occurs at higher T 
than in the CP conserving case (Fig. [7]). Indeed, CP violating effects tend to create an asymmetry in the sterile sector. 
This would push the active system earlier to equilibrium in order to conserve the total null neutrino asymmetry. Since 
L is equilibrated at higher temperature with respect to the CP conserving case, sterile neutrinos are produced more 
efficiently. 

E. e„s = 

In the recent literature, models in which sterile neutrinos mix only with the (mostly) electron ones have been 
discussed as well Q- Hence we also consider a (2 + 1) case in which the mixing angle 6^ s — 0, while 6 es is given 
by Eq. (J6)). The evolution of the sterile neutrino abundance p ss in function of T is shown in Fig. [9l For the sake of 
the brevity, we only consider equal initial neutrino asymmetries L = L e = < 0. We represent the cases L = 
(solid curve), L = — 10 -4 (dashed curve), L = — 10~ 3 (dotted curve), and L = — 10~ 2 (dash-dotted curve). From the 
comparison with the analogous (2+1) case with two active-sterile mixing angles (Fig. 2]), we see that the evolution 
of p ss at different L is qualitatively similar. However, for non-zero asymmetries the production of sterile neutrinos is 
slightly suppressed with respect to the two mixing scenario. 

As in the previous cases, in Fig. [TU]we plot the evolution of the asymmetries Ap a for the different values of initial 
L. In the case L = 0, a value L e = —L s — few x 10 -5 can be dynamically generated. Since u^'s are not mixed with 
the sterile states, their asymmetry remains identically zero. The positive L e can generate a v e -v s resonance in the 
neutrino sector at T ~ 30 MeV. Then, L e becomes negative reaching a value L e ~ — 10~ 4 . After that, the electron 
and the sterile neutrinos go towards flavor equilibrium (with L = 0), reaching it at T ~ 10 MeV. In the case with 
initial L = — 10~ 4 the production of v s starts at T ~ 10 MeV when a v e -v s resonance occurs. We note that since 
i7 asy x A e is a rapidly oscillating function taking both positive and negative values, resonances affect both neutrino 
and antineutrino channels. Then, L e reaches a value ~2x 10~ 3 , while ~ 0. A second series of resonances occurs 
at T ~ 4 MeV, leading L e to zero and L s — — 10~ 4 . In the last two cases, only one resonance occurs at T ~ 3 MeV 
for L = — 10~ 3 , and T ~ 0.5 MeV for L = — 10~ 2 . Once more the sterile neutrino production is triggered by these 
resonances. 



VI. SEMI-ANALYTICAL ESTIMATE OF THE EFFECTS ON BBN 

In order to compute in detail the effects of adding a fourth, sterile neutrino onto BBN, full momentum-dependent 
calculations are necessary. This is essentially due to the fact that the v e and V e distributions enter the weak rates 
regulating the neutron-proton equilibrium and, eventually, the amount of surviving neutrons which will mostly end 
up bound in 4 He nuclei (plus traces of some other light elements), see [8ll. l82j for reviews. We can however provide 
a crude estimate based on a simple physical argument, which has already been used before in this context (see for 
example [48|). Modifying the neutrino sector alters both the Hubble expansion rate and the overall magnitude of 
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FIG. 7: (2+1) case with L = L c — —L^ and i^cp = 0. Left panels: Evolution Ap a = p aa — p aa for the v e (solid curve), 
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the isospin-changing weak rates rj SO , the latter through a change of the neutrino and antineutrino number density 



parameterized here by p ee 3 . Hence the freeze-out temperature Tp, as defined by the condition 

r iso (p ee ,T F ) = H{N cS ,T F ), (37) 

is altered with respect to its standard value Tp — 0.8 MeV due to a higher-than-standard N e g and a lower-than- 
standard p ee . Both effects go in the direction of increasing Tf and, as a consequence, anticipate the freeze-out of 
n/p. This ratio is lower then unity due to the fact that neutrons are heavier than protons by Q — 1.293 MeV, a non- 
negligible quantity compared to the energies involved when T drops to the MeV scale. Since the 4 He mass abundance 
Y p is proportional to the n/p ratio at freeze-out 

Y p oc (-) oc e- Q ' Tp , (38) 



we obtain the estimate 



Y p -T F T F ~ L6 Tp ' [3J) 



3 For the present considerations we assume p ee = ~p ee and neglect the further effect of unbalancing n — > p vs. p — > n rates due to 
asymmetries L, as well as the modified contribution to the Hubble rate due to the asymmetries. For all cases considered here they 
produce only sub-leading changes compared to those illustrated in this section. 
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which confirms that we expect an increase in the produced yield. A simple estimate [48[ for the Hubble parameter 
H oc ^22/7 + N cS T 2 and the weak rates r iso oc (1 + p ee ) T 5 predicts 

1/3 



/ yg2/7 + -V. 



T F a v / ; - , (40) 

V 1 + Pee I 

which immediately illustrates why Y p is comparatively much more sensitive to alteration of the weak rates than to 
the expansion rate via N e g. A perturbative expansion around the fiducial values N e ff = 3 and p ee = 1 yields 

= 0.027 SN eS - 0.17 6 Pee =>■ — £ = 0.044 6N cS - 0.27 <5p ee . (41) 

For specific cases suggested by our previous analysis, e.g. for the (3+1) results of Sec. IV, one finds for example 
that for L = 0, Sp ee ~ while <W e g ~ 1, hence one deduces a variation in the Helium content of 4.4%, which is a 
large number and barely allowed (see e.g. [Hj]). For the largest asymmetries we considered, \L\ = 10~ 2 , the variation 
in JVeff is negligible while Sp ee ~ —0.05, implying again a few percent effects on 1^,. For intermediate values like 
\L\ = 10~ 3 , one expects again effects above the 1% level, this time with both terms contributing. Note that such 
effects are larger than theoretical uncertainties and comparable to observational ones, hence they do imply that sterile 
neutrinos cannot be "easily masked" to BBN: they do have an impact that must be accounted for in any realistic 
analysis combining cosmological observables. The above estimates should be considered only as illustrative, lacking 
a proper account of momentum-dependent effects in weak rates. We plan to provide a more reliable estimate of their 
impact in a future work. 



VII. CONCLUSIONS 



Light sub-eV sterile neutrinos, suggested to solve different anomalies in short-baseline, reactor and solar experi- 
ments [6r]2i i could play an interesting cosmological role providing the amount of extra-radiation indicated by different 
cosmological observations [TtJ ■ However, for the scenarios proposed to fit the different laboratory data, sterile neu- 
trinos would be produced too copiously in the early universe by the mixing with the active species. This would 
create a tension between the laboratory hints and the cosmological observations [24| . A possibility to reconcile sterile 
neutrinos with cosmology is the introduction of a primordial neutrino asymmetry |4lj . that is expected to suppress 
the sterile-active mixing when its strength dominates over the other interaction terms. In this context, we calculated 
the sterile neutrino abundance in the early universe in (3+1) and (2+1) schemes solving the neutrino kinetic equations 
for different initial asymmetries. Considering approximately equilibrium distributions for the active neutrino species, 
the flavor dynamics of active neutrinos plays a sub-leading role in determining the final abundance of the sterile 
species. Therefore, (2+1) schemes are a good proxy for the complete (3+1) situation. Starting with initial neutrino 
asymmetries equal for the two active species, \L\ ~ 10~ 3 would be required to have a significant suppression of the 
neutrino abundance. Otherwise, sterile neutrinos would be produced by resonances between the vacuum term and 
the evolving (oscillating) active v — v asymmetry potential. Opposite initial neutrino asymmetries (hence a globally 
vanishing lepton number) would cuase an enhancement in the sterile neutrino production compared to the above case, 
implying \L\ > 10~ 2 to substantially inhibit their creation. Moreover, in this last case the presence of CP violating 
effects would further increase the sterile neutrino abundance, requiring an even larger initial asymmetry to prevent 
their growth. Both the assumptions of non-dynamical asymmetries and, to a minor extent, of mixing with a single 
active neutrino tend to underestimate the value of \L\ needed for inhibiting the sterile v production. 

Coming to phenomenological consequences, this suggests that some proposed ways to reconcile "hints for a large 
iV" e ff" from CMB with more stringent requirements from BBN, via the introduction of large chemical potentials 
(see e.g. [24j ) are dynamically hard — if not impossible — to achieve. Whenever CMB feels a large N e g due to sterile 
neutrinos of the kind suggested by laboratory anomalies, BBN should feel the same. However, the opposite situation is 
not necessarily true. Even more interestingly, we found that whenever a significant suppression of the sterile neutrino 
production takes place thanks to initial asymmetries, the active neutrinos have partially or mostly decoupled. This 
implies that the (small) fraction of them which oscillates into sterile states is not repopulatcd. Hence, one expects 
different possible regimes: For too small asymmetries, \L\ <C 10~ 3 , the sterile neutrinos are fully populated and their 
"parent" active neutrino spectra are repopulated in the thermal plasma. This implies N e g ~ 4 (for (3+1) scenarios 
considered here) and a tension with cosmological mass bounds, which counteracts the modest fit improvements due 
to a larger N c g. Increasing the asymmetry (\L\ > 10 -3 ) the effect on N e g becomes less and less prominent, and 
completely negligible when \L\ > 10 -2 . However, the lack of repopulation of electron neutrinos would in general 
produce distorted distributions, which can anticipate the n/p freeze-out and hence increase the 4 He yield, to which 
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BBN is much more sensitive than CMB. Finally, for a too large \L\, no production/depletion takes place, but these 
asymmetries in the active flavors would then become an interesting cosmological observable to be associated with 
sterile neutrinos. 

To go beyond semi-analytical estimates, especially to detail the intermediate regime, one has to relax the average 
momentum approximation used in this exploratory study. Due to the momentum-dependence of the resonant con- 
versions between active and sterile neutrinos, a detailed treatment solving the full momentum-dependent equations 
is necessary to derive quantitative phenomenological predictions. We plan to perform this exploration in a forthcom- 
ing article. Also note that our study suggests that the dynamics of sterile neutrinos in the early universe is quite 
dependent on the details of the scenario considered. Pinning down the parameters favored by interpretations of the 
laboratory data in terms of sterile states is crucial in order to treat in the most accurate way only scenarios which 
are phenomenologically attractive. Large scans of parameter space obtained in too crude approximations might miss 
essential aspects of the problem, which is intrinsically non-linear. 
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